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ABSTRACT 

We report the discovery of a statistically significant radial alignment of cluster galaxies in a sample 
of 85 X-ray selected clusters observed in the Sloan Digital Sky Survey. The tendency for alignment 
is a robust feature of our sample and does not vary significantly with individual cluster or galaxy 
properties. We use dynamical arguments to show that a significant fraction of cluster galaxies should 
be undergoing a parametric tidal resonance that can cause their long axes to orient themselves towards 
the center of the cluster potential, and therefore tentatively ascribe the observed radial alignment to 
this dynamical effect. 

Subject headings: galaxies: elliptical — Galaxy: evolution — galaxies: kinematics and dynamics 
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1. INTRODUCTION 

The search for large-scale systematic galaxy alignments 
and their physical causes goes back many years. Even 
before sensitive large-scale photographic surveys existed, 
Brown (1938) suggested that galaxy orientations are not 
isotropic. Hawley and Peebles (1975) were probably the 
first to analyze a large modern galaxy database, finding 
statistically weak evidence that Coma cluster galaxies are 
preferentially oriented toward the cluster center. To our 
knowledge, Coma is the only system where any evidence 
for radial alignment of cluster galaxies has been found 
(Thompson 1976; Djorgovski 1983). On larger scales, 
galaxy alignments have been seen with respect to the 
cluster major-axis and with respect to neighboring clus- 
ters (e.g. Bingelli 1982; Plionis et al. 2003). Galaxy 
cluster data generally support the view that clusters and 
their brightest galaxies do tend to align (c.f. Plionis et 
al. 2003, Struble 1990) although alignment of a signifi- 
cant population of cluster galaxies with the cluster axis 
or toward the cluster center has not yet been established 
(i.e. Fong, Stevenson, Shanks 1990; Trevese, Cirimele, 
and Flin 1992). 

Much work has been done towards formulating plau- 
sible physical mechanisms that could cause alignments 
in cluster galaxies and between clusters. Hawley and 
Peebles (1975) reviewed some of the early speculation 
about how fossil eddy turbulence or primeval magnetic 
fields might be a cause for such anisotropy. Binney and 
Silk (1979) suggested that tidal interactions between pro- 
tostructures could cause anisotropics in galaxy orienta- 
tions, and West (1994) suggests that anisotropic mergers 
along large-scale filaments could induce an alignment of 
the BCGs with the cluster axis and neighbouring large- 
scale structure. Catelan & Porciani (2001) describe how 
tidal interactions can generate a correlation between the 
direction of the angular-momentum vector of a subhalo 
and the surrounding tidal shear field. It is important 
to note that all these mechanisms are primordial, and 
one would expect the alignment they induce at formation 
to decay with galaxy-cluster and galaxy-galaxy interac- 
tions. A particularly interesting question then emerges: 
is galaxy alignment a feature of the initial conditions of 
cluster formation, as assumed in virtually all of the exist- 
ing literature, or is it the result of a dynamical effect that 



grows with galaxy and cluster evolution? These differ- 
ent scenarios lead to distinct alignment tendencies (i.e. 
radial vs cluster-axis vs cluster-cluster alignments) and 
rcdshift evolution. An extensive survey of the strength 
of galaxy alignments in clusters and its correlation with 
individual properties of the parent systems such as mor- 
phology and redshift should help to determine which (if 
any) is the dominant mechanism. 

This letter reports the first detection of a statistically 
significant radial alignment of cluster galaxies for a ho- 
mogeneously selected sample of clusters. In Section 2 
we define our sample and describe the data extraction 
and membership selection procedures. Section 3 presents 
our findings for the resulting samples of galaxies and ad- 
dresses possible systematic errors in our measurement. 
In Section 4 we describe the dynamical tidal interaction 
between a galaxy and the gravitational potential of its 
host cluster and how it can lead to a statistical tendency 
for radial galaxy alignment. The final section comprises 
a brief discussion of these results and a look into the 
future of galaxy alignment studies. 

2. A LARGE SAMPLE OF CLUSTER GALAXIES 

2.1. Sample Definition 

Our data are selected using the largest X-ray and op- 
tical cluster samples currently available in the northern 
hemisphere. We obtain our cluster targets from the Ex- 
tended Brightest Cluster Sample (eBCS, Ebeling et al., 
1998, 2000), a 90% complete, X-ray flux-limited sam- 
ple of the brightest clusters in the northern hemisphere, 
compiled from the ROSAT All-Sky Survey (RASS) data. 
By limiting our sample to known X-ray luminous clusters 
we avoid spurious cluster identifications from projection 
effects and ensure that we are looking at fairly massive 
systems. We cross-correlate the eBCS list with the Sloan 
Digital Sky Survey (SDSS) "footprint" from data release 
3 (Abazajian et al. 2005) and find that, out of 301 eBCS 
clusters, 108 have been observed by the SDSS so far. 

2.2. Data Extraction 

For the data extraction, we define a preliminary cir- 
cular cluster region, 2 Mpc in radius, centered on the 
X-ray position, and a background annulus, 1 Mpc wide, 
4-5 Mpc away from the same point. We assume that this 
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Fig. 1. — Sample galaxy cluster color magnitude diagram with 
fitted red sequence as described in section 2.3 



is far enough for the background region to be unaffected 
by the presence of the cluster but, for most clusters, 
close enough for it to suffer the same sort of system- 
atic uncertainties intrinsic to the survey, e.g. Galactic 
extinction. We assume throughout a flat, A-dominated 
universe, with f2o = 0.3, J7a = 0.7 and h — 0.75. Photo- 
metric and, where available, spectroscopic data are ex- 
tracted for all galaxies that fall within these regions and 
for which position angles have been determined in the 
SDSS. 

2.3. Cluster Membership 

The most practical and accurate way to determine 
which galaxies are gravitationally bound to each cluster 
is to get 3D information directly from spectroscopically 
determined redshifts. Any object with a velocity within 
3cr (where a is the cluster velocity dispersion) of the sys- 
temic velocity of the cluster that is also within 2 Mpc 
projected distance of its center is selected as a cluster 
member (sample A). We use a robust statistical estima- 
tor (ROSTAT) to derive the cluster redshift and velocity 
dispersion (Beers et al., 1990). Even though limited in 
size, this spectroscopic galaxy sample is generated by an 
almost error-free cluster membership selection process. 

In order to increase our sample size we alternatively 
resort to a less accurate, photometric, membership se- 
lection (sample B). Color-magnitude diagrams (CMDs) 
for galaxies in cluster fields allow cluster membership to 
be assessed based on galaxy colors. Figure Q] shows the 
CMD for one of the clusters analysed in this paper (Abell 
1773). The dashed line indicates the best fit to the red 
sequence, a prominent feature of the galaxy population 
of massive clusters. The red sequence is a consequence 
of the -D4000 break - a characteristic of old stellar pop- 
ulations ubiquitous in elliptical galaxies - and Dressier 's 
morphology-density relation (1980). By choosing filters 
which straddle this break at the redshift of the cluster 
(for this work, g' and r'), elliptical galaxies in the cluster 
can easily be selected. 

We use an iterative, automatic fitting procedure for 
delimiting the red-sequence boundaries in order to mini- 
mize observer bias. We collapse all data points along the 
best linear fit to the red sequence and model the resulting 
dispersion histogram as a second-order polynomial (field 
galaxies) with a superimposed Gaussian (red-sequence 
galaxies), as shown in Figure ^ We iterate this process 



until the fit converges, and select cluster members as all 
galaxies falling within ±2<r of the red sequence. In order 
not to miss the cD galaxies, which tend to be bluer than 
the rest, we also extract any bright objects that are bluer 
than the red sequence and within 250 kpc of the cluster 
center. The signal-to-noise ratio for the sample is op- 
timized by taking the extraction radius for each cluster 
to be the point where the local surface density of red- 
sequence galaxies in the cluster area equals that in the 
background region. 

A lower X-ray luminosity cutoff for the host clusters is 
imposed at 0.4 x 10 44 ergs s _1 (0.1 - 2.4 keV) in order 
to exclude systems that are too poor to exhibit a pro- 
nounced red sequence in the available data, reducing the 
final sample to 85 clusters, with redshifts ranging from 
z = 0.02 to z — 0.23. Finally, a control sample (sam- 
ple C) is compiled for each cluster from all the galaxies 
within the background region that fall inside the cluster's 
red-sequence boundaries in color-magnitude space. 

3. ALIGNMENT RESULTS 
Following Struble & Peebles (1985), we describe our re- 
sults in terms of an alignment parameter 8 = ^-—45°, 
where <f>i are the galaxy orientations with respect to the 
cluster center (0° being the radial direction and 90° the 
tangential one) and N is the number of galaxies in the 
sample. For an isotropic distribution, 8 should be con- 
sistent with zero, whereas radial alignment will lead to 
a negative value of 5. We determine 5 for each individ- 
ual cluster and for samples A, B and C and compute 
its lcr statistical uncertainty (derived for an isotropic 
distribution), as = ^== . For galaxy sample A (2210 

spectroscopically selected galaxies in 48 clusters) we find 
8 = —2.21° ± 0.55°, i.e. a weak, but statistically signif- 
icant, preference for radial galaxy alignment at the 4cr 
confidence level. Galaxy sample B (10,472 galaxies from 
85 different clusters) yields only a marginal net alignment 
of S = —0.48° ± 0.25°. However, when limited to galax- 
ies with r' < 18 (the approximate limiting magnitude of 
sample A), this sample shows a much more significant 
galaxy alignment, characterized by 8 = —1.06° ± 0.37° 
(4,853 galaxies). The control sample C (7,396 galax- 
ies over 85 clusters) yields, as expected, a null result: 
8 = 0.69° ± 0.48° and 8 = 0.30° ± 0.30°, with and with- 
out the limitation to r' < 18 galaxies. 

We test the robustness of these results by compar- 
ing the distribution of galaxy orientations in all three 
samples with an isotropic distribution (Fig. and, 
again, find strong evidence for a significant radial align- 
ment of cluster galaxies. A one-sided non-parametric 
Kolmogorov-Smirnov test yields probabilities of 1.7 x 
10" 4 , 3.4 x 10" 3 and 0.35, for samples A, B and C, re- 
spectively, being drawn from an isotropic distribution. 

In order to verify that our measured radial alignment is 
not caused exclusively by a preferred galaxy orientation 
along the cluster axis combined with a flattened galaxy 
distribution along the same direction, we measured the 
global ellipticity, e, of each cluster by computing the spa- 
tial moments of its galaxy distribution. We find no ev- 
idence that the alignment signal is dominated by flat- 
tened clusters: for round clusters (e < 0.25) we find 
8 = -2.48° ± 0.77°, and for flattened clusters (e > 0.25) 
we obtain 8 = -1.92° ± 0.82°. 
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Fig. 2. — Galaxy alignment distribution for samples A,B and 
C. Error bars assume poisson statistics and dotted lines show the 
expected distribution of a 3-dimensional isotropic sample 
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Fig. 3. — The measured alignment parameters 8 and their errors, 
cr for sample A are plotted versus cluster redshift. The dashed line 
shows the sample alignment and its one-sigma error and the dotted 
line shows the expected alignment for an isotropic distribution 



As mentioned in the introduction, any correlations 
of the strength of this alignment with individual clus- 
ter properties will be crucial towards understanding the 
physical mechanisms behind it. To look for a correla- 
tion with morphology, we used a combination of X-ray 
and optical images to separate clusters into 3 morphology 
bins, ranging from extremely relaxed systems to bimodal, 
heavily substructured clusters. We found no statistically 
significant correlation between alignment and morphol- 
ogy or redshift, a result that is perhaps not unexpected 
given the statistical errors in our sample. To see this 
we plot all of our sample A cluster alignment parame- 
ters versus redshift in Figure |3] The dashed line shows 
the best-estimate sample alignment parameter and its 
la error. We find that the reduced chi-squared statistic, 
= 1.10, describes a cluster sample that is character- 
ized by a single value of 6 and which is consistent with our 
previous error estimates. We conclude that further inves- 
tigation into possible correlations of the alignment effect 
with individual cluster properties will require a larger 
cluster sample and/or deeper optical data. 



3.1. Systematics 



We are not aware of a systematic effect that can yield a 
ficticious tendency for radial alignment in samples A and 
B, without also affecting the control sample C. Neverthe- 
less, we delve briefly here into the most obvious sources 
of error. Instrumentally, these are not pointed observa- 
tions, but part of a wide-field survey. The position of the 
cluster centers on the observed fields is therefore random, 
and we can ignore issues such as PSF anisotropics and 
array scanning errors. 

Errors in the SDSS position angles for faint galaxies de- 
termine the limiting magnitude we can work with. The 
SDSS does not provide an error estimate for this quan- 
tity, but we find that the variance in the isophotal ori- 
entation increases steeply around r' = 19.2, and adopt 
this limiting magnitude as an absolute cut-off for the 
red-sequence selection in samples B and C. As discussed 
in section 3., including faint galaxies 18 < r' < 19.2 
in sample B reduces the significance of our result. We 
looked for a dependence on galaxy brightness by split- 
ting the galaxies for each sample A cluster evenly into 
bright and faint subsamples. We found the alignment 
tendency to be effectively equal {^bright = —2.10° ±0.78° 
and Sf a int = —2.37° ± 0.78°), and conclude that for 
r' < 18 there is no evidence that the alignment effect 
is dependent on galaxy brightness. 

4. PARAMETRIC TIDAL EFFECTS ON CLUSTER 
GALAXIES 

In this section we will argue for the existence and im- 
portance of dynamical tidal effects for cluster galaxies 
by drawing an analogy between a cluster-galaxy system 
and the Milky Way (MW)-Dwarf Spheroidal (dSph) sys- 
tems. The dynamical properties of elliptical galaxies in 
clusters are, in some ways, similar to dSphs in the Local 
Group. Fleck and Kuhn (2003-henceforth FK) argued 
that many of the Local Group dSph galaxies have been 
strongly affected by the time- variable MW tide. When 
the internal dSph gravitational timescale is comparable 
to the dSph-MW orbit period important resonance ef- 
fects are likely. The tidal force parametrically excites 
the dSph, and the Mathieu equation describing this inter- 
action exhibits secular growing solutions under a broad 
range of orbit and galaxy density conditions. This is a 
key feature of parametric oscillations, where resonance 
depends strongly on the amplitude of the driving force, 
and weakly on the commensurability of the internal and 
external timescales. FK found that, under dSph con- 
ditions, parametric amplification can have a significant 
influence on the dynamical state of these systems over 
a Hubble time, provided the MW orbital angular fre- 
quency oj — \J GA ^3 1 ~' > ' is within a factor of about 3 of the 
fundamental gravitational frequency of the stellar system 
ljo = TT\/Gp, where M(r) is the enclosed MW mass at 
radius r. This condition is satisfied by some of the local 
dSph, and FK argued that it generally accounts for their 
ellipticities, velocity dispersions, and extended extratidal 
stellar populations. 

Although our SDSS sample spans a wide range of ellip- 
tical galaxy densities and cluster masses, it is worth not- 
ing that a typical elliptical galaxy from this sample and 
its host cluster tend to satisfy the parametric resonance 
conditions. Padmanabhan et al. (2004) used SDSS data 
to generate mass models for elliptical galaxies. Their 
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"typical" SDSS mass model (their Figure 2.) leads to 
a pulsational frequency of luq = 2tt/T with T=30 Myr. 
Cluster mass models based on X-ray and galaxy veloc- 
ity dispersion constraints (e.g. Waxman and Miralda- 
Escude 1995) imply cluster masses of about 1O 14 M at 
a central distance of lOOkpc. This leads to a circular or- 
bit period of 75 Myr which is within a factor of 2.5 of 
the galaxy internal timescalc. The range in masses for 
elliptical galaxies in these SDSS clusters is typically two 
orders of magnitude, and the range in cluster masses is 
similarly large. Thus, our sample includes clusters and 
member galaxies with a range of pulsational and orbital 
periods which broadly encompass parametric resonance 
conditions. 

The importance of the FK mechanism is that it in- 
creases the likelihood of tidal interaction between the 
cluster gravitational field and its constituent galaxies. 
This can have three important effects: 1) it can create 
elliptical galaxies within the cluster, 2) it can generate in- 
tracluster light by stripping stars from the galaxies, and 
3) it tends to align the affected galaxies with respect to 
the cluster center. Equation (25) in FK describes how 
galaxies in nearly circular orbits tend to have stars tidally 
accelerated away from their cores along a line that makes 
an angle 8 with the direction to the cluster center. Here 
tand w — (1 + 2(l^,-oj) )' w ^ n 6 being a a dimensionless 
parameter that expresses the strength of the cluster grav- 
itational tide. Galaxies affected by this mechanism will 
tend to be elliptical and oriented along this direction. 
This equation for 9 shows that, as a galaxy orbits in the 
cluster potential, it will be elongated along an angle that 
is greater (less) than 45° to the cluster center direction, 
depending on whether w is greater (less) than u. Over 
time the effect of the parametric resonance is to decrease 
the galaxy density, so that if luq and to were at any time 
of comparable magnitude then luo would evolve toward 
smaller values as stars are dynamically extracted from 
the galaxy core. While there will also be a tendency 
to circularize the galaxy's orbit as energy goes toward 
heating the stellar system, the change in the orbital w 
is smaller than the density effect on uj and it eventu- 
ally becomes less than to even if it was initially larger. 
This means that those galaxies affected by the FK effect 
in clusters that are at least several crossing times old 
will tend to be aligned along the radial direction (having 
9 < 45°). FK also showed numerically that non-circular 
orbits exhibit the same behavior, as they lead to rotating 
elliptical galaxies with a time averaged preference for the 
galaxy to be radially aligned (cf. FK - Fig. 13). 



It is interesting that a numerical model of cluster galax- 
ies (Muccione and Ciotti 2004) has found that under 
some conditions 10% or more of the total galaxy mass 
can be extracted into the intracluster environment over 
a Hubble time through what they call "collisionless stel- 
lar evaporation" . Because they were simulating galaxies 
with stellar crossing times comparable to cluster orbital 
periods, it is likely that they are also seeing this aspect 
of the parametric tidal mechanism. 

We lack sufficiently detailed cluster galaxy orbit and 
central density information to predict how any given el- 
liptical galaxy should be aligned. On the other hand we 
do expect a statistical tendency toward radial alignment 
of cluster galaxies. This is reminiscent of the dissipativc 
dynamic tidal locking mechanisms that cause planetary 
spin-orbit coupling in solar system objects. 

5. DISCUSSION 

Our analysis suggests that galaxies show a statistically 
robust tendency toward radial alignment in low-redshift 
clusters (z < 0.23). While there may be variations in 
the strength of the effect between different clusters, to 
our measurement sensitivity, the observed radial align- 
ment appears to be a uniform and robust feature of the 
85 X-ray selected clusters we studied using the SDSS 
dataset. The detection of such a net alignment of galax- 
ies in a low-redshift cluster sample seems to point to the 
existence of a mechanism, such as the dynamical tidal 
interaction described above, that can increase the net 
alignment over time so that it remains observable in the 
nearby universe. However, it is also possible that what 
we are seeing are remnants of a primordial " imprint" that 
has not yet been washed away by the tumultuous events 
in a cluster's history. It is vital to expand this analy- 
sis to cover clusters at higher rcdshifts and systems with 
different formation histories so that the effect's evolution 
with time can be studied and the two scenarios can be 
distinguished. With a more sensitive study of the corre- 
lations of galaxy alignment with projected orbit, galaxy 
type, and cluster morphology, we look forward to using 
alignment as a tool to probe, and perhaps distinguish be- 
tween, possible mechanisms behind elliptical galaxy and 
cluster formation and evolution. 
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